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q ; Abstract 

We compute the complete set of Higgs-mediated chargino-neutralino two-loop contributions to 
the electric dipole moments (EDMs) of the electron and neutron in the minimal supersymmetric 
standard model (MSSM). We study the dependence of these contributions on the parameters that 
govern CP-violation in the MSSM gauge-gaugino-Higgs-Higgsino sector. We find that contributions 
mediated by the exchange of WH^ and ZA° pairs, where and ^4° are the charged and CP-odd 
Higgs scalars, respectively, are comparable to or dominate over those mediated by the exchange of 
• neutral gauge bosons and CP-even Higgs scalars. We also emphasize that the result of this complete 

set of diagrams is essential for the full quantitative study of a number of phenomenological issues, 



Or 



P-i, such as electric dipole moment searches and their implications for electroweak baryogenesis. 
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I. INTRODUCTION 



The search for CP-violation (CPV) beyond that of the Standard Model (SM) lies at the 
forefront of nuclear and particle physics. Perhaps the most powerful probes for new CPV 
are searches for permanent electric dipole moments (EDMs) of the electron, neutron, and 
neutral atoms. Null results obtained from these searches have placed stringent constraints 
on CPV in the strong sector of the SM, while present and expected future sensitivities lie 
several orders of magnitude away from expectations based on CPV associated with the phase 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. Various scenarios for CPV connected 
to new physics at or above the electroweak scale naturally imply the existence of non- 
vanishing EDMs that could be observed in future experiments. Thus, a comprehensive 
program of EDM searches could uncover either CPV associated with the "#-term" of the 
QCD Lagrangian, new electroweak scale physics, or both. Each possibility has potentially 
significant consequences for cosmology. The Peccei-Quinn mechanism proposed to explain 
the vanishingly small value of 9 implies the existence of an axion that could account for 
the cold dark matter (CDM), while new electroweak scale CPV could help in explaining the 
observed abundance of baryonic matter through the mechanism of electroweak baryogenesis 
(EWB). 

Among the most theoretically attractive possibilities for new physics is supersymmetry 
(SUSY). SUSY provides an appealing solution to the naturalness problem of the SM. How- 
ever, SUSY has to be softly broken to be consistent with experimental observations. In order 
to solve the naturalness problem, the SUSY breaking scale should be not much higher than 
a few TeV. While the exact mechanism of soft SUSY breaking is not yet known, its effect 
is encoded into the soft terms in the low-energy realization of this scenario. In the minimal 
super symmetric standard model (MSSM), the presence of soft terms implies the existence 
of 40 additional CPV phases beyond the single phase of the CKM matrix in the SM. As 
there exists no known a priori reason for these phases to be suppressed, one expects rather 
sizable EDMs to be generated by one-loop graphs when supersymmetric particle masses are 
below ~ 1 TeV. 

However, the current experiment bounds on electron, neutron EDM, and 199 Hg atom are 
already tight: \d e \ < 1.6 x l(T 27 e cm (90% C.L.) [1], \d n \ < 2.9 x 10~ 27 e cm (90% C.L.) [2], 
and |(i J 4( 199 Hg)| < 2.1 x 10~ 27 e cm (95% C.L.) [3] (For recent reviews of EDM searches and 
their implications for SUSY, see, e.g. Refs. [4, 5]). These results imply CPV phases of order 
10~ 3 or smaller, leading to the so-called "SUSY CP problem". Its resolution, as well as 
that of the related "SUSY flavor problem" , requires some other mechanism for suppressing 
one-loop EDMs (and one-loop flavor changing neutral currents). One possibility is to take 
the masses of the first and second generation sfermions to be of order 10 TeV [6]. In such 
circumstances, the one-loop contributions to EDMs are highly suppressed, and the two-loop 
contributions to EDM, with CP violation from either chargino-neutralino sector or the third 
generation of squarks, may give competitive and even dominate contributions to EDMs of 
the electron and neutron 1 . 

Previous work has considered a subset of these two-loop contributions, including those 
involving third generation squarks [7, 8] and charginos [9-12] whose CPV interactions with 
the gauge-Higgs sector of the MSSM induce an EDM (or chromo-EDM) of an elementary, 



The EDMs of diamagnetic atoms such as 199 Hg will be suppressed in this limit, as they are generated 
primarily by the one-loop chromo-EDM operators. 
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first generation SM fermion. In particular, implications for CP violation at high-energy 
colliders and dominant higher-loop contributions were discussed in detail in Ref. [10]. The 
CP violation from chargino (x + )-neutralino (%°) sector can be propagated to the SM fermion 
though purely gauge boson exchanges. In this case, it has been shown that no CP violation 
can be propagated though 77, 7Z, and ZZ exchanges [11], leaving the WW exchange as 
the only possibility. This contribution was recently calculated in Ref. [11, 12]. CP violation 
can also be propagated through the exchange of gauge and Higgs boson pairs, including 7/1°, 
7 #°, Zh°, ZH°, 7A , ZA°, and WH±. Here, h° and H° denote the neutral, CP-even Higgs 
scalars of the MSSM, with h° being the lightest, "SM-like" scalar; A is the neutral CP-odd 
scalar; and H ± denotes the charged scalars. The contributions due to 7/1°, 'jH , Zh°, ZH°, 
and 7 A exchanges have been studied [9-12]. 

In what follows, we compute the remaining two-loop contributions that survive in the limit 
of large sfermion masses: Barr-Zee [13] type amplitudes wherein chargino-neutralino loops 
communicate CPV to the fermion via the exchange of a ZA° or WH ± pair. We also compute 
the 7/1 , 7#°, Zh°, ZH°, and 7 A contribut ions, and compare our results with the previous 
computations reported in Refs. [9-12]. We report agreement with all previous results. We 
find that, in general, the new contributions are comparable in magnitude to those previously 
computed, or they are even dominant. The ZA° contribution is accidentally suppressed in 
the case of the electron EDM by the 1 — 4sin 2 9w factor, but it is important for the neutron 
EDM. Unlike the case of two-loop diagrams with CP violation from squarks, where it has 
been noted that the 7/1° and 7 A contributions dominate [7, 8], we find - after completing a 
numerical study of the analytic results - that the WH ± contribution is among the dominant 
ones for the electron EDM, and both the Z plus H°, A and the WH^ contributions are 
the dominant ones for the neutron EDM, proving that the inclusion of these contributions 
is indispensable. 

Apart from the implications for EDM phenomenology, our results also have interesting 
consequences for the viability of supersymmetric EWB. Indeed, part of our original motiva- 
tion for computing the loops containing the A and is that the masses of these scalars 
affects the dynamics of supersymmetric EWB. In particular, during a first order electroweak 
phase transition that proceeds via bubble nucleation, the rate at which the neutral Higgs 
vacuum expectation values (vevs) change across the bubble walls depends on m A o (which 
also sets the scale for m H ±). Since the CP- violating asymmetries needed for baryon number 
production are generated during the phase transition by scattering from these vevs, knowl- 
edge of the bubble wall profiles and their dependence on the other MSSM parameters is 
essential for determining the viability of supersymmetric EWB. In general, SUSY EWB is 
enhanced for relatively light m^o - a region in which the corresponding ZA°, ZH° 7A , '-/H 
and WH ± EDM contributions are also enhanced. In our numerical study of the two-loop 
EDMs, we investigate the corresponding m^o-dependence with an eye to these implications 
for EWB. 

Our discussion of these points is organized in the remainder of the paper as follows. In 
Section II we provide details of our two-loop computation and the analytic expressions for 
the results. Section III gives our numerical analysis. We summarize our results in Section IV, 
while additional technical details are provided in the Appendix. We note that, during the 
course of completing our study, a parallel computation of the two-loop EDMs in SUSY using 
an effective field theory approach also appeared [14]. We comment on the similarities and 
differences we find with that analysis. 
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FIG. 1: All the two-loop diagrams with chargino-neutralino loop mediated by Higgs bosons. (Mirror 
graphs are not displayed.) 

II. TWO LOOP EDMS: HIGGS-GAUGE MEDIATED X + -X° CONTRIBUTIONS 

Representative diagrams from the various topologies we consider are shown in Fig. 1. 
For simplicity we have not shown the different crossed graphs or those in which scalar and 
vector boson lines are interchanged. In addition, diagrams involving a photon insertion on 
the SM fermion line are also not shown, as these contributions vanish. Of the remaining, 
non-vanishing diagrams, the contributions with neutral Higgs boson exchange, namely, 7/t , 
7-ff , Zh°, ZH°, and 7 A , ZA°, only involve diagram (a), while the contributions with 
charged Higgs boson exchange W^H^ involve all the diagrams (a,b,c,d). To simplify the 
computation of the latter, we follow Ref. [8] and use the nonlinear gauge [15]. The 
corresponding gauge-fixing term in the Lagrangian is obtained by replacing the ordinary 
derivative that appears in the gauge 

<?. = ~ - idMw^W (1) 

with the U(1)em covariant derivative _D M = <9 M — ieA^ in nonlinear R^ gauge 

C™?^ R t = -±\D,Wt-i£M w cf>+\> . (2) 

Thus, just as the R^ gauge is designed to eliminate mixing between the would-be Goldstone 
boson and the implied by the Higgs kinetic term the nonlinear R^ 

gauge is constructed in such a way that the coupling G ± W T, y arising from the same kinetic 
term is also canceled. As a result, the iy + VF _ 7 vertex is modified from its standard form 
in the conventional renormalizable gauges. 

A direct - and simplifying - consequence of employing the nonlinear R^ gauge is that 
the contribution from diagram (c) vanishes due to the absence of the G ± W 7¥r y coupling. 
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Moreover, an additional simplification can be obtained when carrying out the computation 
in the Landau gauge (£ — > 0). In doing so, one must take care to first compute the £ 7^ 
contributions to the W + W~^ vertex and H^-boson propagators in Fig. 1(b) and carry out 
the appropriate contractions that appear in the one-loop amplitude before taking the £ — > 
limit, since the additional contribution to the W + W~ r ) vertex arising in the nonlinear 
gauge is proportional to l/£. The resulting simplification is that the amplitude from Fig. 
1(d) vanishes as well. This is because the chargino-neutralino loop in (d) is proportional to 
the four momentum of W boson, while the propagator of W is transverse in Landau gauge. 

In carrying out the calculation, we first compute out the one-loop sub-graphs correspond- 
ing to the amplitude for 7(5, ji) — > h(q — t) + g(£, u), where h stands for one of Higgs bosons 
(h°, H°, A , if ± ) having momentum q — £; g denotes one of gauge bosons (7, Z, W) having 
momentum l\ and \i and v denote the vector indices associated with the external photon and 
g, respectively. Gauge invariance implies that the amplitude involves a linear combination 
of the (pseudo) tensors 

p ,u = e ^ qJ p (3) 
= tq u - g^t ■ q . 

For the full two-loop graphs involving the exchange of neutral bosons, only contributes 
in the case of CP-even Higgs exchange, while only T^ v contributes for the graphs involving 
the CP-odd Higgs. Both and T^ u contribute to the two-loop WH ± amplitude. For 
our particular gauge choice, we find that P^ u arises from diagram Fig. 1(a) alone, while for 
W^H* exchange, T^ v requires the sum of both Fig. 1(a) and (b) [graph (b) only generates 
a I ■ qg^ u structure]. We will use these features to explain the origin of the overall, relative 
signs between the various contributions below. 

In obtaining our final results for the two-loop contributions, we use the Feynman rules 
and conventions given in Ref. [16] 2 . We have attempted to express our results in a manner 
that makes easy to directly compare with the earlier work of Refs. [9-12]. We find 



d " S = eQ f C l J' 2 E MD§ aa )M x+ T dx— 1— j(0, -P^), (4) 
M w m 2 s ^ v s > aaJ J x(l-x) Jy 'x(l-x) h y > 
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zs _ ea 2 {Tz fL -2s 2 w Q f )c) m f 



dr = 



16V27T 2 c 2 w s^y M w m 2 

X E lm ( G ab D S,ba ~ G ab D S,ba) M xt f ^Z^S, ^"^"f^ ), (5) 

Jo x xyi x) 



a, 6=1 



A o _ eQ f a 2 c f mf ^ f 1 l-2x + 2x 2 r aA o 

df ~ 8V2n 2 s 2 w M wm 2 A0 ^ ImE - M -iJ dX x(l-x) ^'ifT^' (6) 



2 However, our convention for the Higgs scalar mixing angle a differs from that of Ref. [16]. To facilitate 
comparison with results appearing in the literature, we adopt the convention of Ref. [12]. 
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, ZA o ea 2 (T 3fL -2s 2 w Q f )cj m f 
f ~~ 16a/2tt 2 c^s^ M w m 2 A0 

x£,MG« b E«+G L a A)M xt f dx^j(r ZA o, ^ + ' * )nA ° ),(7) 



/ 



2 ff+ 

ecrc" 



771/ 



2 4 



1 — rr 



j(r WH +, z + ) 

1 — X X 



327r 2 s^c w M w m 2 H+ ^ 
'(lm(MX + M R N R *)M x +x 2 + Im(M> a t + M^N R *)M X .{\ - xf ) 
+ (Im(M> a t - M R N R *)M xt x + lm(M R N^* - M^N R *)M x o{\ - x)) 



(8) 



Here, sw = sm.9 w and cw = cosOw- The S in Eq. (4) and (5) denotes h° and H°. The 
symbol / = u,d,e represents the up quark, down quark, and electron, respectively; Qf 
and rrif are the electric charge and mass, respectively, of fermion /, and T 3 f L is the third 
component of the weak isospin of its left-handed component; finally, j(r,r') is the loop 
function defined in Ref. [11] and given in the Appendix. 
The mass ratios in loop functions are r Zh o 



M 2 z jm 2 



M 2 + /m 2 H+ , 

Xa 

bosons, nih ,h° ,a° ,h- 



M 2 

Xa 

and rw+ 



+ / m h o, r aH o 



M 1 



Xa 



/m 2 H o, r aA o 



M 2 /m 2 h0 , r ZH o 



M 2 +/m 2 A0 , t wh a 



M z /m 2 H0 , r ZA o 



M 2 W 



jm 



2 



M 2 /m 2 H+ , with M w>Zl the masses of W and Z gauge 

> and Al.o > 0, the 



the masses of Higgs bosons, and M + 



and M v o 

masses of charginos and neutralinos, respectively. The coefficients c h u ' A ' H+ and matrices 

D% L H0 ,G R ' L , E R > L , M R ' L , N R,L involve various combinations of the chargino and neutralino 
couplings to Higgs and gauge bosons, and are collected explicitly in the Appendix. 

Before proceeding with our numerical study, we make several comments on the analytic 
results. 

(i) The dependence on the CPV phases in the gauge-gaugino-Higgs-Higgsino sector is 
contained in the imaginary parts of the couplings D R aa , etc. but not separated out 
explicitly. As indicated in the Appendix, these phases arise from diagonalizing the 
chargino and neutralino mass matrices in Eq. (A3). In general, the resulting indepen- 
dent phases are Arg(/xMj6*) and Arg(MjM*), where \x is the supersymmetric Higgs- 
Higgsino mass parameter; (i — 1, 2, 3) are the soft gaugino mass parameters; and b 
is the soft Higgs mass parameter. The SU(3)c mass parameter does not enter the diag- 
onalization of chargino-neutralino mass matrices, leaving two remaining phases. The 
analysis of CPV in this sector is often simplified by assuming that Arg(MiM|) = 0, 
leaving one remaining, independent phase denoted M . In our numerical study be- 
low, we will adopt this simplifying assumption and verify numerically that each of the 
two-loop contributions is proportional to sin0 M . We also comment on the impact of 
relaxing this assumption. 

(ii) The coefficients Cj° and cj° given in Eq. (Al) as well as the matrices D^ L etc. given in 
Eq. (A2) depend in general on tan (5 = v u /vd, where the Vk are the vacuum expectation 
values of the two neutral Higgs scalars, and Z R , which further depends on m A o as 
illustrated in Eq. (A6). This introduces an additional dependence on tan/5 and m A o 



6 



beyond the explicit dependence generated by the Yukawa couplings and dependence 
of the loop functions on scalar masses. 

(iii) The overall sign in the expression for the WH^ contribution is opposite to that of the 
other contributions. The origin of this overall sign can be understood by considering 
the combinations of couplings and Lorentz structures entering the two loop amplitudes. 
To illustrate, we consider the loop of Fig. 1(a) that enters each of the contributions. 
To compare the relative signs of the couplings, we define a general set of interactions 
involving charginos, neutralinos, Higgs scalars, and fermions: 

£v xx = X7 M [AlPl + A R P R ] X V» + ■ ■ ■ 
^xx = x[B L P L + B R P R } X <f> + --- 

Cvei = H [C l Pl + C R P R ] iV^ + ■ ■ ■ (9) 
Ccfu = £[D L P L + D R P R }£<f> + --- 

In Table I below we give the corresponding phases of the couplings obtained from the 
Feynman rules of Ref. [16]. Now we consider the structure of the fermion line in the 
loop, which gives the only other source of a phase difference between the different 
contributions. If I is the momentum flowing through the loop (we may neglect the 
external fermion momenta for this discussion), we have 

ZA° ~ -f [2T 3 P L -2Qsm 2 6 w ] #y 5 
= ^ [-2T 3 P L + 2Q sin 2 6 wl5 ] f 

= [- (2T 3 - 4Q sin 2 6 W ) 7s - 2T 3 ] f 

Zh° ~ Y [2T 3 P L - 2Q sin 2 9 W ] f (10) 

= [(2T 3 - 4Qsin 2 ^) - 2T 37s ] f 

WE ± ~ ^P L fP R = 7 m Pl 9= \^ (1 - 75) f 

For the case of ZA° exchange which involves T^ v from the closed chargino loop, we 
require the 75 term from the lower line to obtain the EDM, whereas for Zh° exchange, 
we have the pseudo-tensor P^ v from the closed chargino loop, necessitating that we 
retain the identity matrix term from the fermion line. For the WH ± exchange con- 
tribution, we require both. Table II gives the resulting overall phase for the various 
contributions. 

TABLE I: Phases of couplings in Eq. (9) as they enter the amplitude of Fig. 1(a) 



Vertex 


ZA 


Zh° 


WH± 


Vxx 


—i 


—i 


+i 


4>xx 


+1 


—i 


+i 


<t)U 


+1 


+i 


—i 


Vit 


+i 


+i 


+i 


Overall 


+1 


+ 1 


-1 
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We observe that the expression for the part of the WH exchange graph arising from 
the T^ v tensor should have an opposite, overall phase compared to the corresponding 
term for the ZA exchange diagram. Similarly, the P^ u component of the WH^ loop 
and the Zh° graph will also differ in overall relative phase. Note that there is an 
additional overall phase that arises between the T^ u and P^ v terms when the identity 

livafi _ _cy- n {iv 
£ a aj3 — 75 

is used in the terms generated by P^ v . Thus, the relative sign between the ZA° and 
Zh° graphs are the same, as are the relative phase between the T^ v and P^ v terms in 
the WPl^ contribution. As we discuss below, the resulting overall sign in Eq. (8) is 
compensated by the signs of various matrix element combinations M^N^f etc. that 
enter the sum over chargino-neutralino intermediate states in the tan/3 ^ 1 regime. 
Consequently, the WH ± and other gauge boson-scalar exchange contributions to the 
electron and down-quark EDMs add coherently in the phenomenologically allowed 
regions of MSSM parameter space. 

(iv) Our results for the 7/1°, 7A , and Zh° amplitudes agree with those of Ref. [9-11], 
including the overall phase. On the surface, our result for the Zh° contribution appears 
to be different from the expression given in Ref. [11]. The difference amounts to 
replacing the x~ l in Eq. (5) by [2x(l — x)] _1 to convert the integral to that of Ref. [11]. 
However, after taking into account the symmetry properties of the integrands in both 
expressions, we have verified (both analytically and numerically) that they agree. 

(v) A direct comparison of our analytic results with those obtained in Ref. [14] is not 
straightforward, since the latter employed an effective field theory approach. We 
note, however, that these authors also include a nonzero result for the ZZ-exchange 
contribution that the authors of Ref. [11] argued should vanish. A direct comparison 
of numerical results is also challenging, since only the dependence of the EDMs on the 
CPV phases was given in Ref. [14], whereas in our numerical study below, we explore 
the dependence on mass parameters and tan/3 for fixed values of 

III. NUMERICAL ANALYSIS 

In this section we numerically assess the impact of the additional, two-loop EDM con- 
tributions discussed above. As mentioned in the Introduction, one motivation for our work 
to consider the complete set of two-loop Higgs-mediated chargino-neutralino contributions 
stems from the framework of EWB. It is therefore natural and well- motivated to focus a 

TABLE II: Summary of signs from fermion line and overall result. The final row is obtained by 
multiplying the overall phases from Table I and the sign obtained from the fermion line. 



Graph 


ZA (T» v ) 


Zh° (P^) 


WH ± (rpHVj 


WH ± (P^) 


7-matrix 


75 


1 


75 


1 


sign 


-1 


+1 


-1 


+1 


Overall 


-1 


+1 


+1 


-1 
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portion of our numerical analysis on a supersymmetric setup which is compatible with that 
framework. Before doing so, however, we investigate the relative importance of the various 
contributions and their dependence on MSSM parameters. 

To that end, we define a benchmark parameter set scenario that will serve as a basis 
for comparison, motivated by the EWB framework, and consistent with phenomenological 
and cosmological constraints, as discussed below. We then proceed with a scan over pairs 
of parameters that govern the size of EDMs, keeping the other parameters fixed at their 
benchmark values. To suppress one-loop EDM contributions, we assume all sfermions to 
be decoupled (we set all sfermion soft breaking masses to 10 TeV, and the trilinear scalar 
couplings to zero, for definiteness) . The gluino mass is entirely unimportant for the phe- 
nomenology discussed here, and is set to 1 TeV. The remaining parameters relevant for the 
two-loop EDM contributions are the absolute values of the gaugino soft breaking masses 
M lj2 and of the higgsino mass parameter //, the heavy MSSM Higgs mass scale (for definite- 
ness we employ here as a free parameter m A o), and tan/3. Our reference benchmark setup 
is defined as follows: 

M x = 145 GeV, M 2 = 290 GeV, fu = 300 GeV, m A o = 300 GeV, tan/3 = 10. (11) 

We consider here one single CP-violating phase, as discussed above. We set this phase 
M = 7r/2, giving the largest CP- violating effect. We verified numerically that the EDMs 
considered here scale proportionally to sin0^ to within an accuracy of 1% . This means that 
(1) all the results we show below can be simply re-scaled when assuming a non- maximal 
CP-violating phase and (2) we show the largest possible size for the EDM contributions we 
consider here 3 . Notice that the values for m h o and the mixing angle a (at the two loop level) 
as well as all the supersymmetric masses and mixing are obtained numerically through the 
FeynHiggs package [18]. 

We have chosen this particular benchmark setup for several reasons. First, this choice is 
potentially compatible with successful EWB. Second, the lightest neutralino relic abundance 
is close to the observed cold dark matter density [19]. If the relic neutralino abundance were 
larger than the cold dark matter density, a mechanism would be needed to dilute the relic 
abundance, with implications for EWB as well [20]. Third, the parameter values given above 
are consistent with collider searches, with precision electroweak data, including the muon 
anomalous magnetic moment, and with the inclusive branching ratio b — > 57 [21]. The 
latter constraint is particularly critical at low m A o, since contributions to b — > 57 from the 
top-quark- loop can be sizable. 

Starting from this reference point, we first illustrate the relative magnitude of two- loop 
contributions in Fig. 2. Here, we explore these contributions to the EDM of the electron 
(left panel) and the neutron (right panel) as a function of tan (3, with all other parameters 
set as in Eq. (11) . To illustrate out results on a logarithmic scale, we show absolute values 
and indicate with black lines positive contributions and with red lines negative ones (the 
only negative contribution we find is from the WW loop, for which we take the expression 
given in Ref. [11]. 

For the computation of the neutron EDM from the quark EDM, we employ the naive 
constituent quark model (CQM) relation 

4 1 

d n = ^d d - -d u . (12) 



3 We remind the reader that EWB implies sin^ M > 10 2 [17]. 
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FIG. 2: A break-up of the various Higgs-mediated chargino-neutralino two- loop contributions to 
the electron (left) and neutron (right) electric dipole moment. Black lines correspond to positive 
values, red lines to negative values. In the right panel we show the various down-quark contributions 
times a factor 4/3, as well as the global up-quark contribution (times a factor -1/3). All SUSY 
parameters, except tan (5, are fixed to the reference setup. 



The computation of d n is subject to considerable theoretical uncertainty associated with the 
non-perturbative strong interaction (for a discussion see, e.g., Refs. [4, 27]. QCD contribu- 
tions to the renormalization group evolution of the quark EDM operators from the weak 
scale to the hadronic scale lead to an overall enhancement factor of 1.53 that multiplies the 
combination of up- and down-quark EDMs on the right side or Eq. (12) [28]. Alternately, 
the use of QCD sum rule techniques to compute neutron matrix elements of the quark EDM 
operators leads to a multiplicative enhancement factor of 1.4 ± 0.6 associated with the qq 
condensate [4]. Ellis and Flores have computed this matrix element by relating the quark 
EDM contributions to the corresponding quark contributions to the nucleon spin as implied 
by polarized deep inelastic scattering measurements, neutron /3-decay, and the Bjorken sum 
rule [29]. This approach leads to a different weighting of the up- and down-quark contribu- 
tions than appears in the Eq. (12) as well as a substantial strange quark contribution that 
is absent from the CQM and QCD sum rule computations. In light of these variations, we 
consider Eq. (12) to provide a reasonable benchmark, bearing in mind that a first principles 
(lattice) QCD computation may yield a different dependence on the light quark EDMs. 
The resulting curves in Fig. 2 lead to several observations. 

(i) As indicated earlier, the sum over all intermediate chargino-neutralino states in Eq. (8) 
compensates for the overall relative sign in the expression for the WH ± contribution. 
After analyzing the individual contributions in detail, we find that although the con- 
tribution from the lightest x + ~X° P air is negative (corresponding to the explicit sign in 
front of the expression), the sum is dominated by the remaining sets of intermediate 
states, many of which carry an opposite relative sign due to the chargino-neutralino 



10 



mixing matrix elements. As a result, the WH ± and other gauge boson-scalar exchange 
contributions to the electron and down-quark EDMs carry the same relative sign. 

(ii) All contributions involving heavy Higgses scale linearly with tan f3, while the WW and 
the contributions involving h° scale as 1/ tan/3. 

(iii) The dominant contributions to the electron EDM appear to be the 7A and 
WH ± loops. However notice that the 7/1 contribution dominates at small tan/5 ~ 2, 
and that the WW contribution is also sizable in that regime. The Z plus Higgs con- 
tributions are suppressed by the T^ eL — 2s'^ v Q e factor, and are relatively subdominant. 

(iv) A similar picture applies to the case of the neutron EDM. Here we explicitly show 
only the break-up for the down-quark EDM contribution (times a factor 4/3), and the 
overall up-quark contribution (times a factor -1/3). As opposed to the electron EDM, 
the ZH°, ZA°, and WH^ 1 contributions dominate the down-quark EDM. 

Using the foregoing considerations, we now study the dependence of the total two-loop 
EDMs as a function of various MSSM parameters. First, we investigate the (tan(3,m A o) 
sector. This sets the mass scale for all loops involving H°, A and H ^ , as well as various 
couplings, directly [see e.g. Eq. (Al)] or indirectly, e.g., through electroweak symmetry 
breaking effects in the neutralino and chargino mass and mixing matrices. We explore 
the EDM dependence on (tan^-m^o) in Fig. 3, setting again all other supersymmetric 
parameters to the values indicated in (11). We indicate in the figure the values of the CP 
violating phase sin</> M such that the resulting EDM equals the current experimental limit, 
d e = 1.6 x 10~ 27 e cm (left) and 2.9 x 10~ 26 e cm (right). For each value of | sinc/> M |, parameter 
space points below the corresponding line are excluded, while those above the line are 
allowed, as indicated. As expected, we find a suppression of the EDM with increasing m A o at 
fixed values of tan /3, thereby allowing for consistency between experimental limits and larger 
values of | sin 0^1 . The overall behavior of the electron and the neutron EDM is remarkably 
similar. Presumably, this similarity indicates that (1) the up-quark contribution to d n is 
subdominant (see Fig. 2, right panel; we actually find that it features significant cancellations 
among the various contributions), and that (2) there are a few dominant contributions to d e 
and to dd that are simply proportional to each other, and hence they contribute in a similar 
fashion. 

To make the connection with EWB, we now analyze the (/x, M 1)2 )-dependence of the two- 
loop EDMs. A generic expectation of the EWB scenario for the MSSM particle spectrum 
includes a relatively light mass scale for the heavy MSSM Higgs sector. This scenario 
depends on the suppression of the net baryon number density generated at the EW phase 
transition with m^o, as pointed out and quantified e.g. in Ref. [22]. It was also realized in 
several analyses (see e.g. Ref. [17, 20, 23] and references therein), that the requirement of 
sufficiently large CP-violating sources for successful EWB prefers a resonant enhancement 
in the higgsino-gaugino sources. This resonance occurs for M\ ~ fi, the resonant neutralino 
baryogenesis funnel, or M 2 ~ /i, the resonant chargino baryogenesis funnel, with M 12 , /i ^ 1 
TeV. An additional resonance could occur for the CPV stop sources, but the latter possibility 
is generally precluded by the LEP limits on the mass of the h°. In both cases, however, 
successful EWB implies sub- TeV masses for higgsinos and gauginos. 

Another generic feature of the MSSM spectrum implied by successful EWB is a light, 
mostly right-handed stop, in order to make the EW phase transition more strongly first 
order. Several studies pointed out, however, that an extended, non-minimal Higgs sector 
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FIG. 3: Exclusion limits for the electron (left) and neutron (right) electric dipole moment, on 
the (tan/3, m^o) plane, for M\ = 145 GeV, M2 = 290 GeV, \i = 300 GeV and for various values 
of the CP violating phase | sm<p^\. For each value of | sin^l, parameter space points below the 
corresponding line are excluded, while those above the line are allowed. 

can also significantly (and perhaps more naturally) enhance the first-order character of the 
EW phase transition [24, 25]. We thus do not regard the requirement of a light stop as a 
necessary feature of a supersymmetric setup giving successful EWB. From the viewpoint of 
EDMs, we considered and evaluated the size of the two-loop stop-mediated contributions 
in [17], and concluded that they are subdominant, provided the left-handed stop is heavy 
enough, as required for EWB in the MSSM 4 . In the interest of singling out the two-loop 
Higgs-mediated chargino-neutralino contributions under investigation here, and in view of 
the above considerations, we do not assume a light right handed stop. 

Having these considerations in mind, we illustrate in Fig. 4 the values of the CPV phase 
giving rise to an electron (left) and neutron (right) EDM equal to the current experimental 
limit, in the (Mi, /i) plane. We assume a minimal supergravity-type relation 5 between the 
gaugino soft breaking masses (where gaugino masses unify at the GUT scale, and their EW- 
scale values are set by renormalization group running), and set here M2 = 2M±. In addition, 
tan/3 and m^o are set to the reference values listed in Eq. (11). In the low \x and low M\ 
region the size of the two-loop contribution exceeds the experimentally viable values (which 
we take to be d e < 1.6 x 10~ 27 e cm and d e < 2.9 x 10~ 26 e cm) for a maximal CPV phase. 
For each value of the CPV phase, we indicate the boundaries of the excluded region. The 
parameter space above the various lines is currently experimentally open even for sin0 M = 1. 



4 Note that this requirement stems primarily from having a Higgs mass consistent with the LEP limits 

5 Notice that the label "[mSUGRA]" obviously doesn't refer to the usual minimal supergravity setup, but 
only to the gaugino mass relation being employed here. The same applies to the label "[mAMSB]" of 
fig. 5. 
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FIG. 4: Exclusion limits for the electron (left) and neutron (right) electric dipole moment, on the 
(Mi,//) plane. We assume the gaugino unification mass relation M2 — 2M\, and vary the CPV 
phase sin^ (see the text for further details on the model assumptions). In the upper right corner 
of the left panel we also show the contour of electron EDM equal to 10 -27 e cm, for maximal CPV 



Electron EDM - sind^=l, M A =300 GeV, M 1= 3M 2 [mAMSB] 



Neutron EDM - sinif =1, M A =300 GeV, M [= 3M 2 [mAMSB] 
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FIG. 5: Exclusion limits for the electron (left) and neutron (right) electric dipole moment, on 
the (M2,/x) plane. We assume here the anomaly-mediated SUSY breaking gaugino mass relation 
Mi ~ 3M2. In the upper right corner of the left panel we also show the contour of electron EDM 
equal to 10 -27 e cm, for maximal CPV phase. 
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In passing, we note that a misalignment of the relative CPV phase between Mi and M2 
does not affect our numerical results. In particular, we find that the main driver for the 
2-loop EDM we consider here is the relative phase between M2 and /i, while a negligible 
contribution originates from the relative phase between M 1 and fj,. This result can potentially 
have profound implications for the interplay between electro-weak baryogenesis and EDM 
searches, which we plan to explore in a future study. 

In summary, our numerical results indicate that the new contributions computed here 
are dominant for the electron and neutron EDM at the two-loop level (in the limit of heavy 
sfermions) and should thus be included in any study of EDMs and CP violation in the 
MSSM. 

IV. SUMMARY AND CONCLUSIONS 

The analysis we have completed gives the EDMs of the electron and neutron in the MSSM 
when the sfermion masses are large, leading to a suppression of the one-loop contributions 
and general dominance of two-loop terms. Considering this regime allows one to circumvent 
the SUSY CP problem associated with the present, stringent EDM limits and sub-TeV scale 
sfermions. It also implies vanishing EDMs for diatomic atoms, assuming they are generated 
primarily by the chromo-EDMs of the quarks 6 . Previous studies of this regime have generally 
also taken all but the lightest, SM-like CP-even scalar to be heavy, thereby suppressing two- 
loop contributions involving the other Higgs scalars as well. In the present study, we have 
not made these assumptions and have, instead, analyzed the dependence of the two-loop 
EDMs on the full gauge-gaugino-Higgs-Higgsino parameter space of the MSSM. 

Our primary result is that contributions arising from exchanges involving one SM gauge 
boson and either the CP-odd neutral scalar, A , or the charged Higgs scalars, H ± , are 
comparable or dominant to previously considered contributions. We have also analyzed the 
prospective implications for MSSM electroweak baryogenesis, whose viability depends in 
part on the values of tan (3 and m A o. We leave a thorough exploration of the MSSM EWB 
to a more comprehensive future study [20]. 
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6 See Ref. [4] for a detailed discussion of the effective operator-dependence of various EDMs. 
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APPENDIX A: COEFFICIENTS, MATRICES, AND LOOP FUNCTION 



The coefficients c„, I, „ ' ' depend on specific types of Higgs bosons and SM fermions 



721 



sin/5' 



d 



cos/3 ' 



sin/5' 



cos/5' 



c u = cot/5, c d = c e = tan/5, c u = cot/5, c d = c e = tan/5. 

R,L 
h°,H°i 



(Al) 



The matrices D% L H0 , G R > L , E R > L , M R > L , N R,L are 
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sin/5(i=Z 2a (^ Sw + Z%c w ) - Z l _ a Z%c w ), 



(A2) 



with tan/5 = v u /vd- 

The Z-t^ are diagonalization matrices of chargino and neutralino mass matrices 
Z T M C Z + = Diag(M x +,M y +), Z T N M N Z N = Diag(M x? , M x o, M x o, M x o), with M y + > 0, 
M o > 0, and 



M, 



a 



( 



( 



M, 



N 



M 2 ^2M W sin/5 
\\f2Mwcosf3 \x 

Mi -M z s w cosp M z s w sinp \ 

M 2 M z c w cosp -M z c w smp 
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(A3) 



The Zr in c h °' H ° and D^ L H0 is the matrix that diagonalize the mass matrix of CP-even 
neutral Higgs bosons 



(A4) 
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which, expressed in terms of a, takes the form 



sum i-os«\ 

cosa smct , 



Z R = I _„„ J • ( A5 ) 



The a, at tree level, can be expressed as 

-(M\ - M 2 z )cos2f3 - y/(M^ + Mjf - AM^M^W 
tana ~ (M% + Mf )sin2/3 ' (Abj 

Notice that a is approximately (3 — n/2 at the limit M A o » M z . Radiative corrections 
modify the latter relation, and we include these effects in our numerical study. 
The loop function j(r,r') is the same as in Ref. [11] 

rLogr j(r) - j(r') 

^ r ) = 73T' -7( r ' r ) = r _ r / • ( A7 ) 
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